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Abstract 



We consider the effects of a massive, unstable neutrino on the evolution of large-scale structure and 
anisotropics in the cosmic microwave background. Comparison with large-scale structure data allows us 
I to rule out a wide range of masses and lifetimes for such neutrinos. We also define a range of masses and 

■^- lifetimes which delay matter-radiation equality and improve the agreement with the data of Cold Dark 

ON Matter models with critical density. 

p PACs Numbers: 98.70Vc, 04.30. +x, 12. 10. Dm, 98.80.Cq 

I 

o 

To appear in Physical Review D. 



1. Introduction 

Large-scale structure is facing a crisis in the post-COBE [1,2,3] era, in that the normalization to cosmic 
microwave background fluctuations on large angular scales is strongly constraining the choice of allowable 
matter power spectra, P(k). In the context of inflationary models, with fi = 1 in mostly Cold Dark Matter 
(CDM) and primordial scale-invariant fluctuations, P(k) is fully specified by a "shape" and normalization. 
The latter is fixed by the amplitude of the temperature fluctuations in the Cosmic Microwave Background 
radiation (CMB) measured by COBE [4]. For models normalized to COBE one finds both the amplitude at 
small scales and the "shape" of the power spectrum predicted by CDM are inconsistent with observations 
of large-scale structure (LSS). 

On small scales the normalization of the power spectrum is usually expressed as as, the ratio of the 
r.m.s. mass fluctuations to the galaxy number fluctuations, both averaged over randomly located spheres 
of radius 8/i -1 Mpc. Observationally the latter has unit amplitude on this scale [5,6], so that as directly 
measures the mass density fluctuations. Observations of pairwise galaxy velocities (on Mpc scales) and 
cluster abundances, effectively probing the power spectrum near 10 Mpc, require as ps 0.5. On slightly 
larger scales bulk flows measured over ~ 10 — 70/i -1 Mpc [7,8,9] would require as fi ' 6 = 1.3(±0.5). Standard 
CDM, normalized to COBE, has a 8 = 1.3 ± 0.1 [4]. 

A phenomenological large-scale power spectral "shape" fitting parameter [10], which basically measures 
the horizon scale at matter-radiation equality, may be defined by (see Eq.(13)) 

where g* counts the relativistic degrees of freedom and g* = 3.36 corresponds to the standard model with 
photons and three massless neutrino species. The large-scale galaxy distribution is phenomenologically fit 
if the "shape parameter" T k, 0.25 ± 0.05 [11], compared with V ~ 0.5 for CDM (with fi = 1 and h = 0.5). 
One can achieve the desired fit in low density models, with either open CDM or ACDM, where g* = 3.36 and 
fio — T/h (such models still have some trouble fitting the normalization of P(k) however [12]). A critical 
density model is also possible if h is decreased [13] or alternatively if g* instead is increased. One means of 
enhancing g*, which we shall explore here, is to fine-tune the mass and decay time of one neutrino species, 
which we shall refer to as v T . 

If the neutrino is massive (1-10 MeV) and decouples non-relativistically, the decay timescale must 
be short (<100s) to avoid excessively perturbing nucleosynthesis [14,15]. (Decaying neutrinos with short 
lifetimes were first proposed as a means of reconciling the normalization on galaxy scales with COBE in 
[16], and have received renewed attention recently [17,18].) For such short lifetimes the early epoch of 
matter (neutrino) domination occurs on such small scales that it is of little interest for seeding any early 
galaxy formation. In contrast, the case of a relativistically decoupling, massive neutrino has no effect on 
nucleosynthesis, and LSS and CMB provide the best constraints. Bardeen, Bond & Efstathiou [19] noted 
that a neutrino with mass in the keV range and a short lifetime (10 4 yr in their example) would provide extra 
power on galaxy scales, and this idea was further pursued in the context of the 17keV neutrino with a lifetime 
of the order of years by Bond & Efstathiou [20]. In this paper we consider the effects of a relativistically 
decoupling v T on the evolution of density perturbations, for a range of masses and lifetimes. Such a neutrino 
changes both the shape of the power spectrum of CDM and the intermediate-scale angular anisotropy in 
the CMB. We have studied in detail the resulting radiation power spectrum for the range of neutrino masses 
and lifetimes which yield a LSS shape parameter in the desired range. We find that decaying neutrinos are 
capable of leaving a distinctive, although subtle, signature in the CMB radiation power spectrum. 

2. Simple Calculation 

The idea that neutrinos may have a mass and be unstable to decay is well motivated in many extensions 
of the standard model of particle physics (see section 5). Indeed if neutrinos with standard annihilation cross 
section have a mass such that 

^2,m v > 92 (tt h 2 )eV , (2) 



they are required to decay so that their energy density today does not overdose the universe. 

In order to build some intuition for the problem, in this section we review and extend a "back-of-the 
envelope" calculation of the matter power spectrum [21,22,23,24]. The interesting feature of a cosmology 
with decaying neutrinos is that, for some range of mass m and lifetime r, the universe has two epochs 
of radiation and matter domination. At very early times the universe is radiation dominated, as in the 
standard cosmology, but for large m there is a redshift, z±, at which the neutrinos become non-relativistic 
and begin to dominate the energy density of the universe. The universe then enters its first matter-dominated 
phase, which lasts until the neutrino decays (into relativistic products) at z dec at which point the universe 
becomes radiation dominated once again. The energy in the decay products redshifts relative to that of 
non-relativistic matter until the universe becomes matter dominated for the second time at z 2 by a non- 
relativistic component, Cold Dark Matter (X) plus baryons (B), with present density £Lnr = &x + &B- 
Assuming the massive neutrinos were relativistic and had standard annihilation cross sections when they 
decoupled from the thermal bath (T ~ lMeV) we estimate 

r iii _1 / 3 

1 + Zl = 2.9 x 10 5 (m/keV) 1 + 0.2(m/keV) 4/3 (r/yr) 2/3 

1 + z dec = 1.3 x lOV/keVr^tr/yr)- 2 / 3 (3) 

1 + z 2 = 2.8 x 10 4 (n NR h 2 ) [l + 0.2(m/keV) 4/3 (r/yr) 2/3 l 

Here we have made the approximations of instantaneous decay at r and sudden change of matter and 
radiation domination and we take the present temperature To = 2.73K. For simplicity we assume the decay 
products stay relativistic until now, thus the bracket in the expression for 1 + z±, which corresponds to the 
ratio of effective entropy degrees of freedom [g* (T±) / g* (To)] is not one (as assumed in [23], for example). 

The constraint on the mass and lifetime from the requirement that the tau neutrino decay products do 
not overdose the universe is 

/ m \ 2 ( t \ „ ^ _ _« ( t c 



< 1.2 x 10 8 ^— — (n h 2 Y , (4) 

VkeW \yij \1.3xlO w yiJ y ' ' KJ 

where whenever necessary we take fio = 1 — ^nr, to = 1-3 x 10 10 yr and h = 0.5. A more restrictive limit 
is obtained from structure formation arguments, which we deal with next. 

The power spectrum we observe today is modified from its "primordial" form by time evolution in the 
expanding universe. The relation between the initial and evolved power spectra is given by the transfer 
function, T(k), (see Eq.(6,7) below). A crude approximation to T(k) may be obtained by realizing that, 
inside the horizon, structure can only grow while the universe is in a matter dominated phase. In this case 
the fractional density contrast Sp/p = 6 oc r/ 2 , with dr\ = dt/a the conformal time. Thus those modes which 
enter the horizon (technically the Jeans scale, but for our purposes the two are equivalent) during a radiation 
dominated epoch have their growth suppressed relative to modes which enter during matter domination. The 
modes which enter the horizon (k = aH) at the epochs of Eq.(3) have comoving wavenumber 

k eq2 = 3.33 x 10- 4 (7iMpc _:L )Vl + 22 

fc dec = 3.33 x 10- 4 (7iMpc- 1 )(l + z dec )/Vl + z 2 (5) 

fc eql = 3.33 x 10- 4 (7iMpc- 1 ) y/(l + Zl )(l + z dec )/(l + z 2 ) . 
Based on this argument we would then estimate 



T(k) 



where 



1 k < k eq2 

(k eq2 /k) 2 fc eq2 < k < k dec 

(Keq2/K de c) K dec < K < fc e qi 

(fc e q2Keql/K d ecK) «eql < K 

P(k) cx kT(k) 2 , (7) 

3 



and the constant of proportionality is fixed by COBE for standard CDM to be about 6 x 10 5 (/i ^^Mpc) . If 
we measure the amplitude of the power spectrum by as, then 



a 8 _ T(k = 0.2 7iMpc _1 ) 
8 (CDM) ~ T(fc = 0.2 7iMpc- 1 ;CDM) 



(8) 



For a limit, suppose we now require that a% not be more than a factor of 10 different from the CDM number. 
This means that 

(k^) 2 (jr")~ 2 °° [ sim P leT «] ( 9 ) 

Here we have assumed that the neutrinos were relativistic when they decoupled, thus this bound applies to 
m < 1 MeV. The bound in Eq.(9) coincides with the usual limit based on structure formation arguments [22], 
obtained by requiring the decay products to become subdominant at recombination, i.e. /9 pro ducts < Pnr at 
t lec . As we shall see, such simple minded power-law approximations to T(k) are not terribly accurate, and 
we shall derive more precise limits below. 

2. Numerical Calculation 

As shown in Fig. 1 approximating the transfer function in the manner outlined above is not quantitatively 
reliable. This is because the concept of a sudden "horizon crossing" is not well justified. To obtain a reliable 
estimate of T(k) requires more than the arguments presented above and we have followed and generalized the 
treatment outlined in [20]. Specifically we have evolved the coupled fluid, Einstein and Boltzmann equations 
for perturbations in a cosmological simulation with decaying neutrinos. The background space time and 
unperturbed energy densities were evolved until the present using 



(10) 



where h = a/a, primes denote differentiation w.r.t. the scale factor a, and dots w.r.t. the conformal time 
T]. We set do = 1 and Ho = 50kms _1 Mpc~ as a compromise between direct measurements and age 
determinations. The evolution is started when all the neutrinos are relativistic so that 

p„ = 2r vPl 
Pt = r„p~ f 

with r„ = (7/8)(4/ll) 4 ' 3 . The pressure of the massive neutrinos is given by 

— = t(x) — — ; — — with x = t (12) 

p Jy ' 3 + xK 1 (x)/K 2 (x) mc 2 V ' 

where the expression with / = 1 is the result assuming Maxwell-Boltzmann statistics (Ki is a modified 
Bessel function) and / is a correction factor fit to a numerical calculation of P/p for Fermi-Dirac statistics. 
The analytic expression gives the asymptotic behaviour for large and small x, making it easy to fit for the 
correction factor. 

For the models in Table 1 we followed the evolution of the perturbations in the synchronous gauge (see 
e.g. [20,25]) treating the baryon-photon fluid as tightly coupled for all time. This causes damping of 6x at 
small scales unless fi# <C 1, but since T(k) is not very sensitive to fi# this limit is justified. Our treatment 
mirrors that of [20] except that we follow the decaying neutrinos from there relativistic to non-relativistic 
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phases. Also the collision term in the Boltzmann equation for the relativistic decay products (which we 
take to be v e and v^) is proportional to A T — 2A„, to be compared with A T — A„ in [20]. The factor 2 
comes from assuming a 2-body final state to conserve momentum, and has almost no effect on the transfer 
functions we compute (the main effect is that it moves the second "plateau" in T(k) to slightly higher k, but 
this "plateau" will not play a role in our studies). When computing the CMB anisotropics (section 4) we 
relax the approximation of tight coupling of the photon-baryon fluid and follow the full photon distribution 
function carefully including the period of recombination. We have compared our results for m = and a 
range of fi# and h with power spectra from Sugiyama [26] and Dodelson [27] and find agreement to better 
than 3% in power (1.5% in temperature). 

3. Large Scale Structure 

We show in Table 1 some measures of large-scale power for 16 models with m and r in an interesting 
range. All models have been normalized to have CMB temperature fluctuations equal to that measured by 
COBE. For the models which are viable, the relative normalization of the large-angle CMB fluctuations 
and the matter power spectrum differ by < 5% from that in the standard CDM scenario. We see that 
for low masses and short lifetimes the results are not very different from normal CDM, whereas for large 
masses and lifetimes they are inconsistent with observations. Using the simple calculations of the preceding 
section we can see that the decaying neutrinos dominated the energy density of the universe at decay only 
if (m/keV) 2 (r/lyr) > 10. Thus for ag, > 1 in Fig. 2 the decaying neutrinos and their decay products never 
dominate the energy density of the universe. 

For the lighter neutrinos no effect is seen until longer lifetimes than we have considered here. We have not 
computed T(k) for models with such long lifetimes, since in this case the models are similar to hot or mixed 
dark matter scenarios where neutrino free streaming causes exponential damping of power on small scales 
[28]. In this case one has to follow the full distribution function for the "light" neutrinos and numerically 
integrate over momentum for each mode, which is a computationally intensive task. Furthermore, because 
of the exponential damping, obtaining sufficient small scale power requires one to have more power at 
~ lOOMpc than LSS would indicate: the LSS data seems to be best fit by a change in the length scale at 
matter-radiation equality (i.e. a change in T), not damping of power on small scales [11]. For the masses 
and lifetimes we have considered the exponential damping sets in at a scales k >few Mpc~ (estimated from 
the horizon size when the neutrinos decay), and does not affect the arguments in this paper. 

We note that for decaying neutrino models with m ~ keV, the second period of matter domination 
leads to more very small scale power than in the standard scenario. This could be useful for early structure 
formation models, however quantitative results are hard to obtain since the details of structure formation in 
the non-linear regime are not well understood. We show in Fig. 3 transfer functions for two of the models 
considered here, the first is a model with T ~ 0.25 and the second is an extreme model from Table 1. Notice 
that only in the extreme case is the plateau from the first period of matter domination clearly visible, and 
even then we find T(k) ~ fc -1 ' 2 rather than constant. 

Efstathiou, Bond & White [12] quote a fitting function 

T(k) = [l + (ak + (bk) 3 ' 2 + (ckfY] ~ 1/V (13) 

with a = (6.4/r)Mpc, b = (3.0/r)Mpc, c = (1.7/r)Mpc and v = 1.13. For a neutrino of mass m w x lOkeV 
and a lifetime r in years they obtain 

"oft (14) 



^0.861 + 3.8(mi O r) 2 / 3 



which is numerically in rough agreement with our results for m ~ lOkeV. Note however that there is not a 
degeneracy between mio and r as this formula suggests, but as shown in Fig. 2,4, it is better to use m 2 T. 
This can be understood as a change in the energy density of radiation after the decay (or g*) which moves the 
second period of matter-radiation equality (see Eq.(3,4)). The extra energy in relativistic products comes 



from the energy "stored" in the v T component while it is nonrelativistic. During this phase (which starts at 
z\ and ends at z^ ec ) the energy density grows compared with that of the relativistic neutrinos oc « ~ t 2 ' 3 . 
Since z\ is determined by T ~ m and a ~ t 1 ' 2 while the universe is radiation dominated (before z\) we find 
the energy is enhanced by an amount depending on (m 2 r) 2 ' 3 . 

Fig. 4 shows the narrow region of masses and lifetimes that give a good fit to LSS, i.e. for which T is 
between 0.2 and 0.3. The models discussed in the next section are within this region. Note that within 
the assumptions we have made these considerations can rule out a large region of (m, r) space which are 
inaccessible to laboratory experiments. In the limit of larger and larger masses, where the lifetimes becomes 
correspondingly shorter, the details of the decay process become irrelevant for all scales of astrophysical 
interest. For the purposes of LSS one has merely increased the number of relativistic degrees of freedom 
in the model and thus delayed matter-radiation equality. One could consider this as an increase in the 
"effective" number of massless neutrino species. We show in Fig. 4 an extrapolation of our results to the 
masses and lifetimes of the rCDM scenario [14]. Note the approximate degeneracy of T with m 2 and r 
appears to hold well over the whole range. 

4. Cosmic Microwave Background Anisotropics 

From the preceding discussion we see that we can fit the LSS data with a combination of masses and 
lifetimes in the keV and year range. The CMB also probes density fluctuations on scales in the 10 — 1000 Mpc 
range, albeit at a much earlier epoch: z ~ 1000. Thus one is led to consider the effect of decaying neutrinos 
on the CMB. In Fig. 5 we show the angular power spectrum of the CMB anisotropics in a model which 
has been adjusted to have the correct shape to fit the LSS data (i.e. T ~ 0.25, [11]). One obvious feature 
is the increase in power on degree scales, which arises due to the smaller sound horizon at last scattering. 
The second thing to note is the observable shift in the secondary and tertiary peaks to higher multipole 
moments (the visibility function in these models is changed only at very low z as expected). This result 
is similar to that of ref. [14] who consider a very massive neutrino which decays on a much shorter time 
scale. In fact, for the models that lie within the band of Fig. 4, we find that explicitly following the decaying 
neutrino, or artificially increasing the number of massless neutrino species to give the same T, lead to almost 
indistinguishable CMB spectra. This is because the range of scales probed by the CMB is large enough that 
all scales were well outside the horizon when the neutrino decayed. 

With the current rate of advance of CMB anisotropy measurements we expect that a measurement of the 
amplitude of the peaks in the CMB spectrum should be possible in the next few years. If the measurements 
continue to prefer a lower peak height then the model in Fig. 5 would be ruled out. To salvage this model a 
lower fi#, tilted initial power spectrum, inclusion of tensor modes or reionization would be required. The shift 
of all the peaks to higher £ is a more robust indication of delayed matter-radiation equality. We have explicitly 
checked that the position of the second and third peaks is fixed by T alone in models with decaying neutrinos, 
artificially increased number of massless neutrinos, and lower h. None of the complications mentioned above 
(changing fi#, etc) would change the position of the second and third peaks, however the amplitude of these 
peaks is dependent on all of them, which makes quantitative comparison with experiment difficult. A further 
technical difficulty is that if fi# is lowered and one wishes to maintain fi^/i 2 at the nucleosynthesis value 
[29,30,31], it becomes necessary to increase h (which would be hard to reconcile with the age estimates) and 
hence find a new m and r which give the same T. 

Decaying neutrino models have slightly more small scale power due to the presence of the first period of 
matter domination, hence it is not implausible to assume that early object formation and reionization could 
have occurred. If this were the case one would expect the anisotropy at small scales to be damped. Note 
that in our "preferred" models however, the amount of extra power on small scales is small, and details of 
structure formation are not well understood. 



5. Particle Physics Models 

Radiative decays, v T — > vj where v stands for a lighter neutrino (y e or v^), in the preferred range 
shown in Fig. 4, are ruled out by several cosmological and astrophysical constraints (see [32,33]). It will 
suffice to mention one of them, based on the absence of photons arriving in conjunction with the neutrinos 
from supernova SN1987A. This limit requires the radiative lifetime of a v T , with mass between 100 eV and 
1 MeV, to be r ra< j > 2.5 x 10 7 s [34]. This still leaves the possibility of invisible decays, e.g. v T — > "iv and 
v T — > v<j), where <f> is a very weakly interacting boson of zero or small mass, most naturally a Goldstone 
boson associated with a spontaneously broken leptonic global symmetry, called Majorons in most models. 

Models with acceptable decay modes of this type were produced [24,35,36] for the heavier v T of [14] 
and for the now dead 17keV neutrino (see for example [37] and references threin). This shows the feasibility 
of models with long lived neutrinos decaying into invisible particles, and in our case model building should 
be much simpler. What considerably complicated the models for the 17keV neutrino was its alleged large 
mixing with the v e , which is not a requirement in our case. In fact the mixing angles that appear in the 
decay rates of the modes just mentioned are in general not the same as those that appear in the mass matrix, 
i.e. are not those measured in oscillation or appearance or disappearance experiments, since operators with 
different flavour structure can contribute. Without this constraint on the mixing from experiment, the v 3 — v e 
mixing in the mass matrix (JJ e z) can be very small (here v 3 is the third mass eigenstate, that we assume 
to be mainly v T ). There is an upper bound on U e3 , coming from neutrinoless double-beta decay, that is is 
not difficult to accomodate. Unless the heavy v T is a Dirac neutrino (and, as explained below, it must be 
Majorana for most of the mass range of interest) it will contribute with amplitude mU 2 3 to the effective 
v e Majorana mass measured in (3(3q v decay. The bound obtained, U 2 3 < (2eV/m), (the (3(3q v lab limit is 
currently m e ff ^ 2 eV) is more restrictive than the present experimental upper bound of about 0.0175 (for 
sin 2 26 < 0.07 and Am 2 > 10eV 2 ; [38]), but not difficult to obtain in models. 

As mentioned above, constraints from energy loss due to the emission of inert right handed neutrinos 
in SN1987A rule out v T Dirac masses in the range O(10keV) to lMeV [39,40,41,42,43]. A related but 
independent bound from SN1987A [44] excludes lifetimes between 3.2 x 10 -17 (m/keV) yr and 1.6(m/keV) yr 
for Dirac masses in the range 1 keV to 300 keV. There is a more dubious bound that would also apply to 
Majorana neutrinos in the same mass range, if the dominant decay is v T — > v e <f>, which rules out lifetimes 
in the range 10 -2 (keV/m) yr to 600(keV/m) yr [37,45]. Energy loss arguments applied to the almost inert 
boson <f>, only require the lifetime to be longer that the duration of the neutrino pulse seen from SN1987A, 
t > 10 (m/ E)s ~ 10 -3 (m/keV)s [35], which becomes relevant only for the heavier v T proposed by [14]. Also 
in contrast to this heavier neutrino, our v T have no effect on nucleosynthesis [15,46,47], and are therefore 
easier to accommodate in particle physics models. 

The decay mode into three neutrinos is more difficult to implement in phenomenologically viable models 
than the decay into a neutrino and a Goldstone boson. For the latter only singlet or mostly-singlet-mixed Ma- 
jorons are allowed by the LEP bound on the number of equivalent neutrino species, essentially N v = 3, mea- 
sured in the invisible decay width of the Z° boson. These models naturally fit well the relatively long lifetimes 
required here. Typical non-minimal singlet Majoron models (and other models, such as familons or com- 
plicated leptonic symmetries) produce lifetimes of order r ~ V 2 /m 3 = 2 x 10 6 (Vyi0 10 GeV) 2 (m/keV) _3 yr, 
where V is the spontaneous symmetry breaking scale (see for example [37] and references therein). The 
desired lifetimes of Fig. 4 correspond to very reasonable values of V around 10 8 — 10 9 GeV. Moreover the 
non-standard annihilations predicted in these models are much smaller than weak (so that the relic density 
of neutrinos is that predicted by the Standard Model, as assumed in this paper). We speak of non-minimal 
singlet Majoron models because in the minimal singlet Majoron model [48], with just one singlet Higgs field, 
the decay rates are suppressed, r ~ V 4 /m 5 [49, 50]. This makes the required spontaneously broken scale 
C(10GeV) leading to phenomenologically more complicated models, if they are viable at all. 

6. Conclusions 

In models with decaying neutrinos the product of neutrino mass 2 and decay time (m 2 r) is strongly 
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constrained by large-scale structure models. The predicted matter and radiation power spectra are approxi- 
mately degenerate in 

— 1/2 

h(l + 0.l[ml v r yT ] 2,3 y (15) 

when flo = 1- We caution however that this formula is only approximate, and it is better to use the 
results of Table 1. Lowering fio or h, or introducing decaying neutrinos, provide a shift in the epoch of 
matter-radiation equality, which improves the agreement of the predicted matter power spectrum with the 
large-scale structure data in CDM models. The latter two scenarios are distinguishable from the former by 
studying the positions of the peaks in the CMB radiation power spectrum on degree scales and smaller. 

In the future we believe that Ho will be measured at cosmological distances via gravitational lensing 
of variable quasars and the Sunyaev-Zel'dovich effect in galaxy clusters, as well as being determined locally 
by a plethora of mostly more indirect measurements. Thus one may be optimistic that, if neutrino decays 
occur in the relevant parameter range to resolve the large-scale structure problem in a CDM model, m 2 T 
can eventually be determined from the CMB. If this turns out to be the case then neutrino physics can be 
studied via temperature fluctuations in the CMB. At present one can already rule out a large range of the 
parameter space from large-scale structure data, see Figs. 2,4. 
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Table 



m (keV) 


t (yr) 


r 


Veo 


v i0 


C8 


0.1 


0.1 


0.51 


304. 


382. 


1.23 


0.1 


1.0 


0.51 


304. 


382. 


1.23 


0.1 


10.0 


0.52 


303. 


381. 


1.22 


0.1 


100.0 


0.52 


300. 


375. 


1.18 


1.0 


0.1 


0.51 


303. 


381. 


1.22 


1.0 


1.0 


0.49 


300. 


375. 


1.17 


1.0 


10.0 


0.43 


284. 


352. 


0.98 


1.0 


100.0 


0.34 


249. 


299. 


0.65 


10.0 


0.1 


0.41 


284. 


352. 


0.97 


10.0 


1.0 


0.29 


249. 


299. 


0.64 


10.0 


10.0 


0.16 


193. 


222. 


0.32 


10.0 


100.0 


0.09 


132. 


145. 


0.13 


100.0 


0.1 


0.16 


193. 


222. 


0.32 


100.0 


1.0 


0.08 


132. 


145. 


0.13 


100.0 


10.0 


0.04 


84. 


89. 


0.05 


100.0 


100.0 


0.02 


54. 


56. 


0.02 



Table 1. Some measures of large scale power for cosmologies with decaying neutrinos, in the limit fi# — > 0. 
We expect ~ 10% variation in these numbers for fi# = 1 — 10%. All models have been normalized to have 
CMB fluctuations consistent with the variance measured by COBE smoothed on 10°. The shape parameter 
T which best describes T(k) is given along with the rms velocity (in km/s) in spheres of radius 60/i -1 Mpc 
and 40/i -1 Mpc after smoothing with a gaussian of 12/i -1 Mpc. In the last column we give the 'normalization' 

<7 8 . 
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Fig. 1: The transfer function (solid) for a standard CDM model (flo = 1) &B = 5%, h = 0.5), along with 
the approximation (dashed) based on Eq.(6). Note that the approximation of a sudden "horizon crossing" is 
not quantitatively accurate. The transfer function is a factor of ~ 3 off at the as scale, indicating the need 
to accurately evolve the perturbations in order to use LSS to set limits on the parameter space. 
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Fig. 2: A contour plot of as in (m,T) space. Contours are for as = 0.5... 1.2 in steps of 0.1, with <r% 
decreasing to the top right of the plot. The lower solid line is <r% = 1.0. We also show a (solid) contour with 
as = 0.1 (upper right). 
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Fig. 3: The transfer functions for a model with m = lOkeV and r = 5 yr (solid) and for m = 100 keV and 
r = 100 yr (dashed). Notice that for the large mass and lifetime the first period of matter domination leads 
to a "plateau" in the transfer function at high k (small scales). For the case m = lOkeV and r = 5yr, which 
has T ~ 0.25 as preferred by LSS, the "plateau" is almost non-existent, though changes in the slope with k 
are visible. 
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Fig. 4: A contour plot of T in (m,T) space which extrapolates to the larger masses and shorter lifetimes of 
the rCDM scenario [14]. The contours are for T = 0.2, 0.25, 0.3 (right to left). Heavy contours mark where 
we have specifically run models, while the lighter contours are the extrapolation to results quoted in [14] 
(the squares). 
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Fig. 5: The angular power spectrum of the CMB. The solid line is for standard CDM (h = 0.5, fio = 1> 
£l B = 5%), the dashed line is for a model with m = lOkeV and r = 5 yr and the same fi and h as sCDM, 
which provides a good fit to the LSS data (T ~ 0.25). Models with p v increased or h lowered to impose 
T = 0.25 are almost identical. Also shown (dotted) is the power spectrum for m = 6 keV and r = 6 yr which 
has T ~ 0.3. Note that the spectrum is quite sensitive to T. 
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